Dynamics of domain walls intersecting black holes 
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Previous studies concerning the interaction of branes and black holes suggested that a small black 
hole intersecting a brane may escape via a mechanism of reconnection. Here we consider this problem 
by studying the interaction of a small black hole and a domain wall composed of a scalar field and 
simulate the evolution of this system when the black hole acquires an initial recoil velocity. We 
test and confirm previous results, however, unlike the cases previously studied, in the more general 
set-up considered here, we are able to follow the evolution of the system also during the separation, 
and completely illustrate how the escape of the black hole takes place. 



Introduction. It is well known that primordial black 
holes and domain walls may have formed in the early 
universe: on the one hand, cosmological density per- 
turbations may have collapsed and formed small black 
holes pj , on the other hand, during the cooling phase 
after the big bang, the series of phase transitions that 
have occurred may have produced extended topological 
structures like domain walls H . In view of the stringent 
constraints that both primordial black holes and domain 
walls provide, understanding how these objects interact 
could help in clarifying many issues regarding the early 
universe and, possibly, be of phenomenological relevance. 

In the past few years, the interaction between black 
holes and domain walls has been object of some study. 
Of particular significance, in relation to our work, are 
the results of Refs. US SHU- Ref. Q considers the 
Dirac-Nambu-Goto approximation for the domain wall 
and proves the existence of a family of static wall solu- 
tions intersecting the black hole event horizon. Subse- 
quent work extended those results to the case of thick 
walls, Refs. HE 011. 

An important problem where the previous studies may 
find application is the scattering of a black hole and a do- 
main wall d . When the domain wall moves towards the 
black hole, it is expected to be captured by the black 
hole and the static configurations mentioned above may 
describe this process in the adiabatic limit, namely when 
the relative motion is very slow. It is likely that the 
membrane will experience large deformations and topol- 
ogy changing processes will take place, but the adiabatic 
limit cannot easily resolve this issue. A full dynamical 
simulation is necessary, and it would also help to clar- 
ify what happens in more extreme situations, which are 
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more likely to have occurred in the early universe, where 
the adiabatic approximation may not be adequate. 

Aside from the motivations mentioned above, the ex- 
citing possibility of observing mini black holes at forth- 
coming collider experiments, as predicted by low scale 
gravity theories 0, 0, 0, 0, 0], has rejuvenated 
interest in the subject. In such models our universe 
has a domain wall structure, the brane, that is embed- 
ded in a higher dimensional space; the standard model 
particles are confined on the brane and gravity propa- 
gates throughout the higher dimensional bulk spacetime. 
Many realizations of this scenario result in lowering the 
Planck scale to a few TeV and share the common pre- 
diction that a small black hole forms when two particles 
collide at sufficiently high energy with small impact pa- 
rameter. Once the black hole is produced, it will emit 
Hawking radiation, and this can in principle be observed 
at the CERN LHC or in cosmic ray facilities. The radia- 
tion of the black hole will partly go into lower dimensional 
fields, localised on the brane, and partly into higher di- 
mensional modes, which will cause the black hole to recoil 
into the extra dimensions. 

The recoil of the black hole in the context described 
above was first studied in Ref. [lfij within a toy model 
consisting of two scalar fields, one describing the black 
hole and the other a possible quanta emitted by the black 
hole in the process of evaporation; the coupling between 
them is fixed as to reproduce the probability of emis- 
sion of a scalar quanta from the black hole as prescribed 
by the Hawking formula, and a delta function potential 
is used to mimic the interaction between the black hole 
and the brane. In the approximation that the interac- 
tion with the brane is negligibly small, it is shown that, 
as soon as a quanta is emitted in the extra dimension, the 
black hole will slide off the brane. Although this has the 
important observable effect of a sudden interruption of 
the radiation on the brane (Ref. explores some pos- 
sible phenomenological consequences), it is not clear how 
the separation process occurs. To clarify this issue and 
to have further evidence for the conclusions of Ref. [la , 
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two of the present authors have considered the interaction 
between 'mini' black holes and branes from the different 
perspective of studying the dynamics of branes in black 
hole spacetimes, Ref.|l7|. Specifically, by treating the 
brane in the Dirac-Nambu-Goto approximation, we have 
shown that, once the black hole acquires an initial re- 
coil velocity perpendicularly to the brane, an instability 
develops and the brane tends to envelop the black hole. 
This suggests a mechanism of escape for the black hole 
due to the reconnection of the brane. In some approxi- 
mation the time of escape can be estimated and it was 
found to be shorter than the evaporation time. This fact 
may have important phenomenological consequences. 

The previous claim certainly deserves reconsideration 
and one of our motivations is to test the robustness of 
the conclusions of Ref. [lTj by modeling the brane as a 
domain wall composed of a scalar field. Also, unlike the 
Dirac-Nambu-Goto case previously studied, modeling the 
brane as a domain wall allows one to describe reconnec- 
tion phenomena. In our case, this means that we can 
follow the evolution of the system while the separation 
of the black hole occurs, and completely illustrate how 
the escape takes place. 

Domain wall dynamics. The system we intend to study 
consists of a domain wall intersecting a black hole. As 
customary, we consider black holes whose mass m ex- 
ceeds the fundamental Planck scale M, in order to make 
quantum gravity corrections negligible. Also, the size 
of the mini black hole is assumed to be smaller than 
the characteristic length of the extra dimensions. We 
shall restrict our analysis to the case when the wall ten- 
sion is small enough, as this allows us to ignore the self- 
gravity of the wall. In this regime, the mini black hole 
is completely immersed in the higher dimensional space- 
time and can be adequately described by asymptotically 
flat solutions [II El] 

ds 2 = -f(r)dt 2 + /- x (r)dr 2 + r 2 d^ 2 , (1) 

where dQd~2 is the line element of a (d— 2)— dimensional 
unit sphere, 

and rjj is the horizon radius. We use units in which 
c = h=l. 

Since the character of vacuum structures is generally 
insensitive to the details of the model, we limit ourselves 
to a scalar effective field theory with quartic potential 

V($) = ±(&-r) 2 ) 2 . (2) 

The domain wall solutions of this model in flat space 
have tension a ~ \f\rf and thickness t ~ 1/yXrj. The 



coupling constant A is of order A 4_d , with A < M be- 
ing the cutoff of the effective field theory. We can ne- 
glect the quantum corrections to the domain wall-type 
configurations at weak coupling A^ d ~ 2 ^ 2 ^§> r\. This 
implies a t 1_d , and hence we have ar'jf 2 /M ^> 
{r H /t) d - 1 {M d /M) { - d -^^ d - y >, where we introduced the 
d-dimensional Planck mass M d = G 1/(d - 2 \ Requirement 
of negligible self-gravity therefore implies that the wall is 
not very thin and/or the mass of the black hole is suffi- 
ciently large. 

The evolution equations describing the dynamics of the 
domain wall intersecting a black hole, are given by 

= -r 1 (r)9 t 2 $ + -^9 r (r d - 2 /Ma,$) 

+ 27 ■ TOi ft [(sme) d ^d e t>] - V'(*) , (3) 
r^{smO) a 6 

where we assumed 0{d— 2)-rotational symmetry around 
the axis perpendicular to the equatorial plane of the black 
hole. 

The initial conditions are specified so as to mimic the 
initial recoil of the black hole. For this purpose, we im- 
pose that <& and 3$ at t = are given by a static kink 
profile in the flat spacetime boosted in the direction of 
the symmetry axis 

x / „x , I M rcos6> — vt \ 

^ v (t,r,9)=r ] t a nh[^-r l —==-\ , (4) 

with v constant. 

Boundary conditions have to be specified on the outer 
boundary, on the symmetry axis and on the horizon. At 
'infinity', since the gravity of the black hole is expected 
to be negligible, we assume that $ reduces to the flat 
boosted form The regularity conditions at the sym- 
metry axis and on the horizon do require 

d e $(t, r, 6) = , at 9 = 0, tt , 
d t ${t,r,0) = , at r = r H . (5) 

In order to solve Eq. JSJ numerically, we use a mixed 
spectral and finite difference method. The first step is 
to decompose the general solution in terms of a complete 
basis of smooth global functions [2(J. Here, for conve- 
nience, we adopt the Chebyshev polynomials of second 
kind, and write the general solution as 

N 

$(t,r,6) = J2<Pn(t>r)U n (cos9) , (6) 

n=0 

where N is the number of harmonics used. The boundary 
condition along the symmetry axis is trivially satisfied, 
as the Chebyshev polynomials U n (cos9) are regular at 
9 = 0, tt. The decomposition JBJ is used as ansatze in 
equation ©, and some trivial algebra allows us to recast 
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Eq. in matrix form 

ip = f 2 (r)<p" + f(r) 



f'(r) 



d-2 



f(r) 



where we have defined 
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(l-.T 2 ) i 



-dx, 



U(x)V'(<P)\/l- x 2 dx 



(7) 

(8) 
(9) 



x = cos 9, and the vectors tp and U are given by 

ip = (^o(^,r'),¥'l(^,? , ), • • • ,<PN(t,r)) , 
U = {U Q {x),Ui{x),--- ,U N (x)) . 

Equation Q is a system of coupled two-dimensional non- 
linear partial differential equations, which can be solved 
by using finite difference methods. 

We define ip~ , ip and <p + , which represent the con- 
figurations at three subsequent time steps, i_ = t — St, 
t and t+ = t + St, respectively. We now define a one- 
dimensional grid along the radial direction, and label the 
points as r^. Then the values of a function g(t,r) on a 
grid point at t = i_ , t and t + will be denoted by g~ , 
Qi and gf. In terms of the coefficients of the expansion, 
the initial conditions are written as 



and 
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$„(0, rj,x)U(x)\/l - x 2 dx 



$ v (5t,n, x)\J{x)\Jl - x 2 dx 



(10) 



(11) 



where $> v (t, r, x) is the boosted domain wall solution, 
Eq. JIJ. As for the boundary conditions, we have 



<p(t,r > Too) 



&v{t, r > Too, a;)U(a;) - \/ 1 — x 2 dx . 



at the outer boundary r^, and 
<p(t,r = m) 



, 



on the horizon. 

To solve the system of Eqs. J7J) we use a central dif- 
ference discretization scheme. In this way the R.H.S. of 
(J7J) can solely be expressed in terms of (p, whereas the 
L.H.S. will depend on <p~ and ip + . Therefore, once the 
boundary and initial conditions are specified, the system 
of equations is completely determined. The poten- 
tial term F is evaluated by first reconstructing the field 
profile $ and then use it to evaluate the function V'(&). 



As last step, V'(Q) is convoluted with the basis functions 
according to formula ©. 

The details of the simulations will be reported else- 
where, here we briefly illustrate the main result. Fig. ^ 
shows six snapshots describing the evolution of the en- 
ergy density of the domain wall, 



p(t, r, 6) 



2/(0 



4> 2 



— ^ 2 e + Vm. (12) 



The plot presents the result of a simulation obtained by 
using N — 30 harmonics and 100 grid points in the ra- 
dial direction. Standard numerical checks have been per- 
formed and the evolution has been carried out in various 
cases with the thickness of the domain wall ranging from 
0.01 to 1 times the size of the horizon radius and for vari- 
ous initial velocities, v ranging from 0.01 to 1. The results 
obtained in the case of a domain wall confirm the sugges- 
tion of Ref. 0, 0] : the domain wall envelops the black 
hole, that completely separates from the mother brane; 
after the pinching occurs, the mother brane relaxes into 
a stable configuration. 

Modeling the brane as a domain wall, we have given 
further evidence in support of the conclusions of Refs. 0, 
llH Iri| that the black hole intersecting a brane may es- 
cape once it acquires an initial recoil velocity. As in the 
previous work, we have neglected the self-gravity of the 
wall, but the present treatment allows us to describe re- 
connection processes of the domain wall. By following 
the evolution of the system during the separation, we suc- 
ceeded in illustrating completely how this process takes 
place. 

There is one important question which has not been 
considered here: what happens when the tension is non- 
negligible? Intuition suggests that the back reaction may 
prevent the separation at least for small initial velocities, 
and an energy barrier has to be overcome to induce in- 
stability leading to pinching. In other words, when the 
tension of the brane is switched on, it might be possi- 
ble to find a sequence of static configurations before we 
encounter a critical unstable configuration above which 
the instability sets in. However this problem is rather 
complicated and beyond the scope of the present article. 
Study to understand this issue is in progress and we hope 
to report on this soon. 
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FIG. 1: Evolution of the domain wall. The figure shows how 
the energy density of the wall evolves, explicitly illustrating 
the separation process and escape of the black hole. The 
reported simulation refers to the five dimensional case, d = 5, 
and a domain wall whose thickness is 1/3 of the Schwarzschild 
radius and v = 0.3. 
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